Statistical properties of the continuum Salerno model 
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The statistical properties of the Salerno model is investigated. In particular, a comparison between 
the coherent and partially coherent wave modes is made for the case of a random phased wave packet. 
It is found that the random phased induced spectral broadening gives rise to damping of instabilities, 
but also a broadening of the instability region in quasi-particle momentum space. The results can 
be of significance for condensation of magnetic moment bosons in deep optical lattices. 
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Nonlinear wave propagation models are of importance 
in a wide variety of settings, such as in nonlinear optics 
p], , in water wave propagation Q , and in plasma sys- 
tems 0. Moreover, also discrete systems show important 
nonlinear properties that give rise to solitary structures 
in e.g. lattices (see, e.g. [5| and references therein for 
a review). The discrete Salerno model |(| combines on- 
site and intersite nonlinear terms in a lattice in order to 
encompass the different characteristics of nonlinear inte- 
grable and non-integrable lattice models. It has also been 
discovered as an asymptotic form of the Gross-Pitaevskii 
equation, with applications to Bose-Einstein condensa- 
tion of atoms with magnetic moments in deep optical 
lattice traps 0, @ • Recently, the continuum limit of the 
Salerno model with focusing onsitc and defocusing inter- 
site nonlinearities was investigated @ and soliton solu- 
tions were found. 

In this Brief Report, we will investigate the modu- 
lational instability properties of the general continuum 
limit of the Salerno model. Furthermore, the effects of 
random phased wave packets on the stability properties 
will be studied using the Wigner formalism. Results re- 
garding the effect of random pulse phases will be given. 

The Salerno model was derived as a quantum modi- 
fied discrete nonlinear Schrodinger equation, giving the 
time evolution of the field amplitude on the lattice. The 
relevant equation takes the form |f| 

J a t $„+($„+i+$„_i)(i+ A1 |$„| 2 )+2^|$„| 2 $ n = 0, (i) 

where $„ = $>(t,x n ) is the field amplitude on the nth 
lattice site located at x n and v and fx are the real param- 
eters determining the strength of the onsite and intersite 
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nonlinearities, respectively. The continuum approxima- 
tion is introduced by the ansatz $(t, x) — *£(t, x) exp(2ii) 
and using the expansion ^ n ±i ~ $ ± d x *ff + d^.^f/2 the 
continuum limit of the Salerno model is obtained from 
Eq. (1) according to 

id t $ + 2{v + /i)|*| 2 * + (1 + /i|*| a )Sg* = 0, (2) 

where VP is the continuum correspondence of the complex 
lattice field amplitude. We see that by letting p — ► 
we obtain the nonlinear Schrodinger equation, where the 
sign of v determines the character of the soliton solu- 
tions 0, while v — > gives a continuum version of the 
Ablowitz-Ladik equation 

Next, we derive a wave-kinetic equation for the wave 
function We first introduce the Wigner function for 
the wave function <3> according to jlOl fllL Il2| 



p(t,x,k) = — 
Ztt 



dZe^{H>*(t,x + Z/2)*(t,x-Z/2)), 

(3) 

where the angular bracket denotes the ensemble average 
The Wigner function is a generalized distribution 
function for the quasi-particles representing 'J, and it 
satisfies the normalization 



I(t,x) = (\V(t,x)\ 2 ) = / dkp(t,x,k) 



(4) 



Applying the time derivative to Eq. J3J and using Eq. 
P)l we obtain the wave-kinetic equation 

dtp + 2kd x p + 2/{ sin (i d x dk) [2(i/ + p) + p(k 2 + dl)} 
+ cos (i dxdk) pkd x }p = 0, (5) 

where the sin and cos operators are defined in terms of 
their respective Taylor expansion, and the arrows denote 
the direction of operation. Equation (J5J describes the 
propagation of partially coherent wave modes taking on- 
site as well as intersite nonlinearities into account. 
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Next, we look for unstable perturbations around the 
homogeneous solution po(k). We let p(t, x, k) = po(k) + 
Pl(k) exp(iKx — idt), where \p\\ <C po, and linearize (4) 
with respect to p\. Equation (JSJ with the intensity Q 
then gives the nonlinear dispersion relation 



1 



dk 



[2(u + p) + pk 2 ][p (k + K/2) ~ Po (k- K/2)] 
- 2K(1 + pl )k 

r (6) 

The dispersion relation 10 generalizes the results of pre- 
vious modulational instability analysis of the nonlinear 
Schrodinger equation, where the Wigner formalism 0] 
as well as the mutual coherence method 0] has been 
used. 

For a coherent background spectrum, i.e. po(k) = 
IoS(k), the dispersion relation JHJ reduces to 



±{(l + /i/ ) [(l + |/xIo) K 4 - 4(v + p)I K 2 ] } 



1/2 

j\ ' ^ ' L \ Z ' "/ v ' / " J| 

(7) 

We see that for the nonlinear Schrodinger case (p = 0), 
we obtain the standard growth rate T = —itt = K[Ai/I — 
K 2 } 1 / 2 , showing the existence of a modulational instabil- 
ity for focusing [y > 0) nonlinearity. 

For a random phase background wave function, we 
may represent the background quasi-particle distribution 
function by the Lorentz distribution 



Po(k) 



it k 2 + A 2 ' 



(8) 



where A is the width of the distribution function. Us- 
ing the Lorentz distribution in the nonlinear dispersion 
relation, we obtain from (6) 



ft = -2iAK (1 + \pl ) 

±[{l + ph) [(1 + f/i/o) K 

+2A 2 K 2 pI Q (1 + i(Ml ) } 



4 - 4(i/ + p)I K 2 ] 

1/2 



(9) 



We note that the dispersion relation reduces to Q 
when A — > 0. Moreover, the case of a pure Kerr nonlin- 
earity (p = 0) gives rise to the a linear damping term for 
the modulational growth rate. In general, the broaden- 
ing of the background quasi-particle distribution is seen 
to give rise to a damping of the modulational instability 
growth rate through the first term in the dispersion rela- 
tion @ , but through the intersite nonlinearity parameter 
p we also have a new coupling due to the finite width of 
the distribution function. 

Next, we numerically investigate the properties of the 
growth rate T = Im(f2), deduced from (8). In Fig. 1 we 
have plotted the modulational instability growth rate for 
different parameter values. We have normalized the pa- 
rameter v to one, and used A = and A — 0.1. For 
p = A = we retrieve the standard result for the modu- 
lational instability growth rate in the case of the nonlin- 
ear Schrodinger equation. This behavior is, however, af- 
fected by a finite width in the background wave spectrum. 
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FIG. 1: The growth rate F plotted as a function of the per- 
turbation wavenumber K for different parameter values. All 
curves have normalized the parameter v to one and intensity 
Io = 0.25. The heavy solid curve has /i = A = 0, and cor- 
responds to the coherent nonlinear Schrodinger modulational 
instability; the thin solid curve has fi = and A = 0.1; the 
heavy dotted curve has p = 0.5 and A = 0, while the thin 
dotted curve has /i = 0.5 and A = 0.1; finally, the heavy 
dashed curve has /j, — —0.5 and A = 0, while the thin dashed 
curve has p = —0.5 and A = 0.1. Thus, all cases show damp- 
ing through the broadening of the background wavenumber 
distribution. 



All cases show a Landau-like damping of the growth rate, 
when a finite width of the background distribution po is 
introduced. Thus, appropriate random phasing of the 
wave function could act as a useful means of stabi- 
lizing, e.g. Bose-Einstein condensates in deep optical 
traps. It is interesting to note that there is an interplay 
between the onsite and intersite nonlincaritics through 
the coefficients v and p, respectively. As expected, the 
case of competing nonlinearities, i.e. pv < (here repre- 
sented by a repulsive intersite nonlinearity), gives rise to 
a reduced growth rate as compared to the case of a Kerr 
type nonlinearity (p = 0), since any perturbation growth 
in one nonlinearity can be decreased by the other nonlin- 
ear term. However, for p, v > (the dotted curves) we 
note a significant increase in the growth rate. 

To summarize, we have performed a wave-kinetic study 
of the Salerno model in the continuum limit. A Vlasov- 
like equation has been obtained for the quasi-particles, 
representing the dynamics of the wave function "J" in 
phase space. A modulational instability analysis has 
been carried out, and a comparison between the coher- 
ent and incoherent cases has been made. It was found 
that the interaction between the onsite and intersite non- 
linearities and the finite width of the background distri- 
bution function gives rise to significant changes in the 
modulational instability growth rate. 
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